Is the 130 GeV Line Real? A Search for Systematics in the Fermi-LAT Data 
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Our recent claims of a Galactic center feature in Fermi-LAT data at approximately 130 GeV have 
prompted an avalanche of papers proposing explanations ranging from dark matter annihilation to 
exotic pulsar winds. Because of the importance of such interpretations for physics and astrophysics, 
a discovery will require not only additional data, but a thorough investigation of possible LAT 
systematics. While we do not have access to the details of each event reconstruction, we do have 
information about each event from the public event lists and spacecraft parameter files. These data 
allow us to search for suspicious trends that could indicate a spurious signal. We consider several 
hypotheses that might make an instrumental artifact more apparent at the Galactic center, and find 
them implausible. We also search for an instrumental signature in the Earth limb photons, which 
provide a smooth reference spectrum for null tests. We find no significant 130 GeV feature in the 
Earth limb sample. However, we do find a marginally significant 130 GeV feature in Earth limb 
photons with a limited range of detector incidence angles. This raises concerns about the 130 GeV 
Galactic center feature, even though we can think of no plausible model of instrumental behavior 
that connects the two. A modest amount of additional limb data would tell us if the limb feature is 
a statistical fluke. If the limb feature persists, it would raise doubts about the Pass 7 processing of 
E > 100 GeV events. At present we find no instrumental systematics that could plausibly explain 
the excess Galactic center emission at 130 GeV. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

The search for non-gravitational signatures from 
WIMP (weakly interacting massive particle) dark mat- 
ter has generally been approached from three different 
directions: missing energy searches at colliders, direct 
searches for the recoil of nuclei from underground detec- 
tors, and indirect methods including searching for dark 
matter signals from cosmic rays (CR) and multiwave- 
length astronomical observations [iHg]. 

For indirect detection, distinguishing the dark mat- 
ter signal from conventional astrophysical backgrounds is 
challenging (for a recent review on indirect searches with 
gamma rays see 0]). Among various possible signatures, 
gamma-ray line emission is a long-sought "smoking gun" 
for dark matter annihilation [f|, as no plausible astro- 
physical background can produce such a line signature. 1 
Gamma-ray line(s) could be produced by dark matter de- 
cays or annihilations into two photons, or two-body final 
states involving one photon plus a Higgs boson, Z bo- 
son, or other neutral non-SM particle. In most models, 
the branching ratio to lines is loop suppressed relative to 
the continuum emission, and one would have expected 
to see the continuum first in e.g. MSSM models [e.g. 



1 A narrow feature is possible in theory [see [^] . 



[Tfjj . Although this theoretical prejudice led most pre- 
vious studies to focus on continuum searches, there are 
models being proposed that allow high line to contin- 
uum ratios [e.g . Q, [TlHl5| . However, previous searches 
in EGRET [Of and Fermi-LAT data (xMll did not find 
any indications for a gamma-ray line signal and presented 
only upper limits on the line flux. 

First indications for a spectral feature around 130 GeV 
were found by Bringmann et al. [20] | in context of virtual 
internal Bremsstrahlung signals from annihilations. The 
first claim for a significant line at the Galactic center 
(GC) was made by Weniger [2l|. Both works focused 
on spectral fitting to photon events in regions of interest 
in the inner Galaxy designed to maximize S/N. Weniger 
found a line structure with A.6a (3.2a after the trials fac- 
tor correction) at 130 GeV, and argued against an obvi- 
ous instrumental cause. This claim was quickly followed 
up and disputed by a number of groups |22l |23| . 

Subsequent work by Su & Finkbeiner approached the 
problem with template fitting, which takes into account 
the spatial distribution of events along with spectral 
information, assuming various profiles (Einasto, NFW, 
Gaussian) for the DM distribution |24|. If the template 
is correct, this allows extraction of the DM signal with 
higher S/N. This work found 6.6a (5.1a after the tri- 
als factor correction) for an Einasto profile centered 1.5° 
west of the Galactic center, and also suggested that there 
may be two lines, at about 111 and 129 GeV. The lower 
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energy line is tantalizing because it matches the expected 
energy of a Z^j line if the higher energy is the 77 line. 
These findings have inspired a number of models and fur- 
ther analysis of the Fermi data [H EE HUM HESi • 

Recent evidence for lines at 111 GeV and 129 GeV with 
a local significance of 3.3cr from Fermi unassociated point 
sources suggests an annihilation signal is present [45| [but 
see H(| , as does the claim of line emission from galaxy 
clusters at 130 GeV [13] ■ Neither of these would stand 
on their own, but they provide support for the hypothesis 
that the Galactic center line signal is produced by dark 
matter annihilation. 

The high statistical significance of the line feature mo- 
tivates a search for systematic errors in the LAT data 
that could mimic a line in the Galactic center. Confir- 
mation by Imaging Air Cherenkov Telescopes like HESS- 
II might be possible as early as next year [48j], but in 
the meantime a thorough study of LAT systematics is 
urgently needed. We do not have access to the details 
of the reconstruction of each photon event, which would 
allow us to study how it developed in the tracker and 
calorimeter. However, we do have information about 
each event from the public event lists and spacecraft pa- 
rameter files. We can use this information to search for 
any line-producing artifacts in the detector frame, and 
investigate if they could map onto the Galactic center. 

The Earth's atmosphere provides a convenient source 
of photons for systematics tests. The continual cosmic- 
ray cascades in the Earth's atmosphere produce gamma 
rays with dN/dE ~ E~ 2 8 [49j |, Because these so-called 
'Earth limb photons' result from atmospheric cascades, 
they are produced by interactions in a highly boosted 
frame, and cannot contain line emission. 

The rest of the paper is organized as follows: In Sec- 
tion fll] we briefly define key parameters of LAT photon 
events and describe the survey strategy. We examine 
peculiarities of the Galactic center observation, possible 
systematics of the LAT in the instrumental frame, and 
study whether they could fake a 130 GeV line signature. 
In Section [TTT1 we concentrate on a suspicious subset of 
Earth limb photons that shows a line-like excess at 130 
GeV. We search for correlations with the GC line events, 
and introduce an energy remapping model as possible ex- 
planation for spurious signals. Finally in Section llVl we 
discuss our findings and what is required to clarify the 
status of the 130 GeV excess. 



II. THE 130 GEV EXCESS 

In this section, we briefly summarize the standard sur- 
vey strategy of Fermi and define the basic parameters of 
each event. We address the question of whether obser- 
vations of the Galactic center are peculiar in a way that 
could enhance instrumental effects towards this direction 
and search for suspicious trends in other regions of the 
sky and the reconstruction parameters. 



A. Standard survey strategy and definitions 

With a field of view of ~ 2.5 sr, the LAT can survey 
the entire sky in two orbits. In standard survey mode, the 
LAT points north of zenith towards the orbital pole by an 
angle Z mt ± on one orbit, and south of zenith by the same 
angle on the next orbit. In this mode, the LAT pointing is 
confined to the plane perpendicular to its orbital velocity. 
The slews are performed with a repeating pattern of 17 
waypoints defining Z Ioc ^ as a function of time. 2 Z m( ± 
was 35° at the start of the nominal mission on August 
4, 2008 until May 7th, 2009, and was changed to 50° 
on September 3rd, 2009 for better thermal management 
of the downward-facing battery radiator. 3 This rocking- 
angle profile, combined with the precession of the orbit 
every ~ 53.4 days, allows the LAT to observe the whole 
sky with approximately uniform coverage. 

Fermi spends over 95% of the mission time in standard 
survey mode. This is only occasionally interrupted for 
pointed observations of targets of opportunity (ToOs). 
During such times the LAT may point at larger zenith 
angle than usual, even at the horizon. Fermi's survey ob- 
servations are halted during passages through the South 
Atlantic Anomaly, resulting in an exposure differential 
between north and south of ~ 15%. In addition, sur- 
vey mode is occasionally interrupted by Autonomous Re- 
points of the observatory for triggered gamma-ray burst 
follow-up observations, and for calibration. 

The reconstructed arrival direction of photons in celes- 
tial coordinates, LAT coordinates, and Earth coordinates 
is described by parameters in Table fl] 8 is the recon- 
structed incidence angle of the photon event with respect 
to the LAT boresight (defined as the +z axis). The +x 
axis is the line normal to the Sun-facing side of the space- 
craft, i.e. the solar panels, which are parallel to the y axis, 
face roughly the +x direction. (f> is the azimuthal angle of 
incidence with respect to the +x axis. The Zenith angle 
Z is the angle between the reconstructed event direction 
and the zenith line, which passes from the Earth center 
through the satellite center. All angles are in units of 
degrees. 

Due to the increased rocking angle Z roc k~ 50° since 
September 2009, photons from the Earth limb entered 
the FOV of the LAT: At a spacecraft altitude a, the 
geometric (unrcfractcd) horizon is seen at zenith angle 

Z hoi = cos" 1 (-J^—) + 90° . (1) 
The Fermi orbit is nearly circular with 535 km < a < 



2 The survey rocking angle profiles are available at 
http://fermi.gsfc.nasa.gov/ssc/observations/types/allsky/ 

The effective dates and times for each profile are provided. 

3 Various profiles with rocking angle 39° , 40° , and 45° respectively 
have been tested for relatively short periods, including a 3 orbit 
profile test that overweights the south. 
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Parameter 




Range 


Description 


name 


min 


max 




e 





~ 80 


Polar coordinate (instrument frame) 


4> 





360 


Azimuthal coordinate wrt +x (instrument frame) 




-1 


1 


Cartesian coordinates; oc (sin 6 cos <f>, sin 6 sin <f>, cos 6) 


Z 





~ 113 


Zenith angle (horizon at Z = 113°) 


I 


-180 


180 


Galactic longitude 


b 


-90 


90 


Galactic latitude 







180 


Angle to Galactic center; cosi/> = cos £ cos b 


^rock 


-110 


110 


Rocking angle (boresight angle N of zenith) 



TABLE I: Event parameter definitions. The Cartesian coordinates (x, y, z) are defined such that the z-axis corresponds to the 
LAT boresight, and the j/-axis is parallel to the solar panels. 



564 km, yielding Z\ lor in the 112.7° to 113.3° range, with 
the tangent point some 2400 km distant. At this distance, 
the ~ 100 km height of the atmosphere subtends about 
2°, or roughly 111° < Z < 113°. Combined with the 
large rocking angle, the Earth limb events are dominantly 
seen near the incidence angle 9 = Z\ lor — Z IO( ± ~ 62°. 

In Table [TT| we define the event samples used through- 
out this work: 'Standard events' are all events with the 
zenith angle cut Z < 100° recommended by the LAT 
team to exclude Earth limb photons (see Section lLTH) ; 'In- 
ner Galactic plane' refers to a part of the Galactic disk 
close to but without the center; 'Galactic center' events 
come from a radius of 3° around £ = b = 0°; 'Earth limb 
events' have zenith angles Z > 110°, and are completely 
dominated by photons generated in CR cascades in the 
atmosphere. 'Line' events refer to subsets with energies 
between 120 and 138 GeV. This energy range is selected 
since the dominant line at the GC is found to be around 
129 GeV [Hi, and the FWHM of the relevant LAT energy 
dispersion is about 13.6% at that energy (2lj . 

Throughout, we will use P7CLEAN_V6 events from Aug 
4th 2008 to September 5th 2012, with 10 GeV < E < 
500 GeV and the good-time- interval cuts DATA_QUAL==1 
and LAT_C0NFIG==1, however without the commonly 
adopted cut on the rocking angle |^ roc k| < 52°, unless 
otherwise stated. This last cut would remove low inci- 
dence angle Earth limb events, which will be of special 
interest below. 

In Fig. [TJ the gray dots show the distribution of all 
events (i.e. without the Z cut) with energies above 100 
GeV as a function of the instrumental incidence angles 8 
and <f). The contribution from the Earth limb is clearly 
visible at 6 > 60°. For comparison, the red and blue dots 
show the 'GC line' and the 'Earth limb line' events that 
will be discussed below. 



B. Peculiarities of the Galactic center observation 

The fact that the dominant 130 GeV line signal is near 
the Galactic center raises a number of concerns. The 
gamma-ray flux at the GC is somewhat brighter and 
might have a harder spectrum than neighboring regions. 



gray: all £>100 GeV events 90° 




270° 



FIG. 1: Incidence angle distribution of all events with E > 
100 GeV (gray dots; see discussion in Sec. Ill A} , of GC line 
events (red dots; see Table [TXJ and of suspicious Earth limb 
events (blue dots; see TableM 'Earth limb line'), as a function 
of the instrumental coordinates 6 and <f>. The Earth limb 
contribution is clearly visible in the gray dots at 8 > 60°. 
Note that just because a GC line event was observed at high 
6 does not mean it was observed near the horizon. 

Also, the GC is near the ecliptic (/? » —5°) and is ob- 
served in a restricted range of angles in instrument coor- 
dinates (9, (f>) when the Sun passes near it. We consider 
whether these facts could exacerbate any systematic er- 
rors in the LAT data to produce a spurious signal. 

1. Hypothesis: The Galactic center is bright, so 
instrumental artifacts are more significant there. 

The hardware trigger rate of the LAT is typically about 
10 3 -10 4 Hz, with the rate of accepted SOURCE and CLEAN 
class events below 3 Hz [50|, and the rate of E > 100 GeV 
events orders of magnitude lower. In light of these low 
trigger rates (and assuming steady sources and CR back- 
ground fluxes) the LAT instrumental response cannot de- 
pend on the brightness of an observed region. Further- 



Sample 


Cuts 


iV(> 100 GeV) 


JV ^ J> ±UU uC V J 

JV(>30 GeV) 


AT 1 '-if HI T" 1 \ 

iV(>100 GeV) 


Standard events 


Z < 100° 


5093 


13.4% 


9.6% 


Inner Galactic plane 


Z < 100°, 3° < \£\ < 30°, |6| < 2° 


703 


16.9% 


9.8% 


Galactic center 


Z < 100°, ip < 3° 


82 


17.4% 


9.8% 


Galactic center line 


Z < 100°, i> < 3°, 120 GeV < E < 138 GeV 


26 






Earth limb 


Z > 110° 


3120 


10.2% 


9.2% 


Earth limb line 


Z > 110°, 30° < 9 < 45°, 120 GeV < E < 138 GeV 


45 







TABLE II: Definition of six samples of events used throughout this work. The number with E > 100 GeV is given for each 
sample, along with the indicated ratios. All samples have 10 < E < 500 GeV. 



Events from different test samples 




FIG. 2: Energy distribution in various event samples as dis- 
cussed in the text. None of them show the excess events 
around 130 GeV seen in the Galactic center. The red dotted 
lines indicate 111 and 129 GeV. 



more, at E > 100 GeV, the Galactic center is only mod- 
estly brighter than the surrounding regions, so that re- 
lated effects should also appear away from the GC. Oth- 
erwise, fake 130 GeV events would have to be mistakenly 
mapped from either lower energy (E^IO GeV) gammas 
or much lower energy photons (e.g. X-rays from the IE 
1740.7-2942 microquasar [U or 511 keV photons (H^) in 
which the Galactic center is much brighter. It is difficult 
to see how this could happen. 

Bright regions provide samples with a high gamma- 
ray-to-CR ratio, which are used for calibration pur- 
poses by the LAT team [H(J. Besides that, their main 
virtue is that they feature a large number of events so 
that the impact of an instrumental effect, like e.g. en- 
ergy reconstruction or acceptance anomalies that affect 
the reconstruction of gamma-ray events, can be statis- 
tically more significant there. We check as a warm-up 
whether we find indications for suspicious features at 
130 GeV. Besides the Galactic center (with an inten- 
sity of ~ 8 x 10" 8 ph cm" 2 s" 1 sr" 1 above 100 GeV), 
we consider the Earth limb (with all incidence angles 
and an intensity of ~ 3 x 10~ 7 ph cm~ 2 s _1 sr _1 above 
100 GeV [49]) and the inner Galactic plane (~ 9 x 
10~ 8 ph cm" 2 s _1 sr -1 ), see Table HU The sample with 
the highest number of celestial photons (but a smaller 
intensity, ~ 5 x 10 -9 ph cm~ 2 s _1 sr _1 , and hence a 
higher CR contamination) is simply the whole sample 
of standard events excluding the region around the GC 



(i/j > 10°). For these samples, we show the correspond- 
ing energy distributions in Fig. [2] None of the sam- 
ples exhibit an 0(1) excess at 130 GeV as observed in 
the Galactic center. Note that in the vicinity of the 
bright point sources Vela and Geminga (within 0.8° ra- 
dius, corresponding to the 95% containment angle) only 
four E > 100 GeV events were measured, at 102.5, 111.5, 
123.8 and 205.5 GeV. 

In Fig. |3] we show fits to the energy spectra of the GC, 
inner Galactic plane, and Earth limb. The model fits 
a line on a power-law background, convolved with the 
instrumental response, as in [2l| . The TS value is defined 
as TS = — 21n£ nu n/£ a it, where £ nu ii( a it) is the likelihood 
of a fit without(with) a line. The line normalization is 
constrained to be non-negative. While the TS value of 
18.8 in the GC is relatively large, there is no indication 
for line contributions in the Inner Galactic plane or Earth 
limb samples. 

We emphasize that signal significances calculated from 
the GC region as defined in this paper do not represent 
the full significance of the putative signal, which is higher 
in regions with optimized signal-to-noise ratio [20l [211 ] or 
when extracted by a template analysis 24j. We use the 
GC region as defined here since it should be dominated 
by line events, making it a good starting point to look 
for suspicious trends in the data. 



2. Hypothesis: The Galactic center has a hard spectrum, 
making energy mapping errors more significant. 

The Galactic center black hole (Sgr A*) is visible up to 
20 TeV, and other sources in the GC may be unusually 
hard. If these high energy photons are occasionally mis- 
reconstructed with energies close to 130 GeV, this could 
produce a line feature in the data that would appear 
preferentially at the Galactic center. In Table [TTl we list 
the ratio of E > 300 GeV events to E > 100 GeV events 
and of E > 100 GeV to E > 30 GeV events for our 
samples. The GC spectrum is not much harder than the 
rest of the Inner Galactic plane, but the latter shows 
no sign of a feature at 130 GeV (Fig. [3]). Note that 
about ~ 20 events contribute to the central part of the 
line feature seen at the GC Assuming that the GC 
spectrum is not harder than dN/dE oc E~ 2 (normalized 
to the number of events above 150 GeV), the GC events 
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GC, TS = 18.8 




E [GeV] 

Inner Galactic plane, TS=0.0 




200 



E [GeV] 

Earth limb, TS=0.2 




3. Hypothesis: the GC observations have a restricted range 
of incidence angles on the instrument. 

If the Galactic center were predominantly observed at 
specific angles in LAT instrumental coordinates, associ- 
ated instrumental problems could be projected onto the 
Galactic center simply for geometric reasons. In late De- 
cember (June), the Sun crosses the Galactic disk, and 
the angular distance between the Sun and the Galac- 
tic center (anti-center) decreases to ~ 5°. Since Fermi 
keeps the solar panels aligned to the Sun (the Sun is in 
the x — z plane), this leads to an increase of Galactic 
center events at instrument azimuth angles of <$> ~ 0° 
(cf> ps 180°). This behaviour is clearly visible in the up- 
per left panel of Fig. 0] where we show the (f> distribution 
of > 10 GeV events from the Galactic center as a func- 
tion of time of year. The distribution in the upper right 
panel reflects the precession pattern with a ~ 53.4 day 
period. Integrated over time, the <fi and 9 distributions 
look like in the lower panels of Fig. 21 where they are 
shown as a function of the Galactic longitude. Close to 
the GC, the <j> distribution becomes significantly bimodal. 
Note that this effect occurs along the full ecliptic, but the 
intersection with the Galactic disk close to the center is 
accompanied by the largest number of events. 

It is tempting to relate the location of the 130 GeV 
excess in the Galactic plane to this inhomogeneous (p 
distribution. As a test, we select standard events from 
the full sky (excluding ip < 10°) in the range <f> = 
-20 ... 20° mod 180°. The total number of these events 
is ~ 10 times larger than the number of Galactic cen- 
ter events; an anomaly in the event reconstruction near 
these </> angles would appear in this data sample with high 
significance. However, the energy distribution shown in 
Fig.[S]shows no significant feature at 111 or 129 GeV. 



E [GeV] 



FIG. 3: Spectral fits to the GC, inner Galactic plane, and 
Earth limb samples. The green line shows the null model (a 
power-law), whereas the red line shows the alternative power- 
law + line fit; the dotted blue lines are the two components of 
the alternative model. The red (black) dotted lines indicate 
129 GeV (13.6% FWHM around 129 GeV). Note that the 
significance found in the GC region does not represent the 
full significance of the putative signal. 



cannot come from energies above 300 GeV, since there 
would not be enough events to make up the 130 GeV 
excess even if all of them were incorrectly mapped to 130 
GeV. The possibility of energy mapping errors of events 
close to 130 GeV will be discussed in Sec. Mil 



C. Peculiarities of the instrument 

Instrumental effects are likely to be correlated with 
specific instrumental coordinates rather than sky coordi- 
nates. In this subsection we search for suspicious trends 
at 130 GeV as a function of the event incidence angles, 
quality parameters, and arrival times. We also search 
for the possibility of 'hotspots' at other regions of the 
sky; they could indicate instrumental effects or an exotic 
source population. 

An error in the LAT effective area (e.g. various cuts 
at ~ 130 GeV are less/more restrictive for some reason) 
could in principle explain the GC line, but would also 
produce a line in the other samples listed in Tab. [TXJ 
Exactly this already happened in the past at energies 
close to 10 GeV However, no 130 GeV feature is 

seen in the main test samples outside of the GC (though 
in a small subsample that we will discuss below). The 
LAT team estimated the point-to-point correlation of the 
effective area on scales relevant for line searches to be on 
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Binning plot of 0<1()° standard events Binning plot of v><10° standard events 




[deg] l[deg] 



FIG. 4: Upper panels: <j> (left panel) and 6 (right panel) distributions of standard events close to the Galactic center, as a 
function of time modulo one year (with Jan 1st at the origin). In Dec (Jun) the Sun passes near the Galactic center (anti- 
center). In these periods, most GC events are observed at (j> f» 0° ((/> « 180°), since the LAT x — z plane is determined by the 
Sun to keep the solar panels oriented. Lower panels: <f> (left panel) and 8 (right panel) distribution of standard events along 
the Galactic disk. Close to the GC, the distribution becomes significantly bimodal. The red and blue dots are as in Fig. Q] 



the level of 2% [53[, which is not enough to explain the 
0(1) GC or Earth limb features. 

If the GC line is due to non-rejected CRs, it is hard 
to understand how mono-energetic particles could be 
present in the CR spectrum. What is more, a bright 
region like the Galactic center has a large gamma-ray-to- 
CR ratio, such that a CR contamination would affect it 
last, not first. 



1. Hypothesis: A 130 GeV features is visible at specific 
incidence angles on the instrument. 

One concern is that a 130 GeV feature in the recon- 
structed events is visible only for events with certain inci- 
dence angles (9, (f>) in instrumental coordinates. The 130 
GeV excess would be then imprinted on regions of the 
sky that are predominantly observed at these problem- 



atic angles. 4 

To study line-like features at different incidence angles, 
we split up the {9, (f>) plane in different regions as shown 
in the left panels of Fig. [6j We then analyse the spectrum 
of all events that hit the detector at these incidence angle 
patches separately and search for lines. The fits are per- 
formed as in Ref . [2lj , but the energy window is fixed to 
80 to 200 GeV, we assume a flat acceptance and for a line 
we take a simple Gaussian with a standard deviation of 
6% (to approximate the central part of the non-Gaussian 
LAT energy dispersion). 

The left panels of Fig. [6] show from top to bottom the 
results obtained (1) using all events, (2) using Earth limb 
events only, and (3) using standard events only; the latter 
two are disjoint subsets of the former one. For each (9, <j>) 



4 A study of the x-y positions of the events inside the tracker and 
calorimeter is not possible using public data only. It is however 
difficult to conceive how these parameters could be correlated 
with certain regions of the sky. 
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Standard events close to x—z plane and ip>\a 




1.8 2.0 2.2 2.4 2.6 



FIG. 5: Energy distribution of standard events away from 
the GC with incidence angles close to the LAT x — z plane 
(perpendicular to the solar panels). The red dotted lines in- 
dicate 111 and 129 GeV; the thin red line shows GC events 
for comparison. 

range we quote three numbers: the significance of a 129 
GeV line-feature in terms of the TS value, the signal- 
to-background ratio of the putative line signal, and the 
number of events above 80 GeV that contribute to this 
region. The TS values are also indicated by color for 
better visibility. 

We also compute the TS as a function of energy in 
each region of (9, <f>) to check whether there is anything 
unusual about 130 GeV (Fig.jBJ right panels). In particu- 
lar the lower two plots which correspond to independent 
data sets show clearly that the enhanced TS values at 
different energies in different panels are not correlated. 
Lastly, we note that the significance of the largest TS 
value, TS ~ 11 (at 115 GeV), has a p-value of about 
~ 0.3 when taking into account ~ 2 x 23 x 2 trials over 
a x1=2 distribution (two Z ranges, 23 incidence angle 
panels, ~ 2 independent search regions [54j|). 

Fig. [6] suggests that the TS values observed in some 
panels in case of Earth limb (Z > 110°) events are well 
within the statistical expectations once the large number 
of associated trials is taken into account. However, it can- 
not be excluded that the 130 GeV excess in these parts 
of the Earth limb indicates a real instrumental effect. 
Then the energy reconstruction at different incidence an- 
gles would have to somehow additionally depend on the 
zenith angle Z of the events (e.g. via a subtle correla- 
tion with the rocking angle). We postpone an in-depth 
discussion of this possibility to section IIIII in any case 
more Earth limb data will help to verify or falsify this 
possibility. 

2. Hypothesis: The Galactic center events are flagged as 
badly reconstructed. 



the event-level analysis, which inform about the qual- 
ity of calorimeter and tracker event reconstruction. 5 
CTBCORE describes the probability that the direction 
estimate is good (roughly the probability that the recon- 
structed direction falls within the nominal 68% contain- 
ment angle), CTBBcstEnergyProb the probability that 
the reconstructed energy falls into the core of the energy 
dispersion [50J . We show these parameters for the the GC 
line events in the upper panel of Fig. [7] the background 
histogram shows the distribution of the these parameters 
in standard events above 100 GeV. No significant bias in 
the distribution appears. 

The red bars in the lower panel of Fig. [7] indicate the 
first tracker layer that shows evidence of a particle hit for 
the best track reconstruction (TkrlFirstLayer) in case of 
the GC line events. Tracker layers are 0-17, where 
is closest to the calorimeter and 6-17 (2-5) corresponds 
to FRONT- (BACK-)converting events (tracker layers and 
1 do not have conversion foils due to requirements of 
the three-in-a-row trigger primitive). The dark gray bars 
show the distribution averaged over all > 100 GeV stan- 
dard events. The light gray bars show for comparison 
the distribution of a 'dirty' > 100 GeV event sample. It 
is defined as all SOURCE-CLEAN events with |6| > 5°, and 
hence has a high CR contamination. The distribution of 
the latter is significantly biased towards the first and last 
tracker layers, whereas the distribution of GC line events 
is compatible with the expectations for standard events. 



3. Hypothesis: There are 'hotspots ' with line-like features 
in other sky regions. 

The 130 GeV excess is located in a narrow range of 
about ~ 5° around the Galactic center. The existence 
of other significant 'hotspots' (localized regions with a 
significant line-like emission at or around 130 GeV) in 
the sky could indicate an unexpected instrumental ef- 
fect, or point to a new exotic source class [23] . In ab- 
sence of knowledge about the distribution of such possible 
sources, we adopt the following strategy. We select par- 
tially overlapping regions of 6° x 6° size, centered along 
the Galactic disk at b = 0°, \£\ = 13.5°, 16.5° . . . 178.5°. 
The size of the region is chosen to include roughly 50-100 
events above 80 GeV (like in our 'GC region'); regions 
close to the GC are excluded to avoid contamination with 
the 130 GeV excess. In each region we calculate the 
TS value for the presence of a line at various energies. 
The details of the fit are the same as above for Fig. [3] 
(in particular the energy window is always fixed to 80- 
210 GeV). The distribution of TS values obtained in this 
way is shown in Fig. [8] by the dotted colored lines for 
different line energies; a combination of TS values for all 
line energies is shown by the solid black line. Its tail 



Although we do not have access to all details of the 
event reconstruction, the extended LAT event files con- 
tain a few figure-of-merit quantities from the first step of 



http: //fermi .gsfc .nasa. gov/ssc/dat a/ analysis/document at ion/Cicerone/ 
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FIG. 6: Left panels: Significance for a ~ 130 GeV line- like excess in different parts of the 6-cj> plane (for 9 — 0-80°). From top 
to bottom, the results are derived from all events, from the Earth limb events, and from standard events. The three numbers 
show the TS value for the presence of a 129 GeV line, the signal-to-background ratio and the number of events above 80 GeV 
inside the considered 9-<f) region. #-cuts are made at 15° , 30° , 45° and 60° . Right panels: The significance for a line-like excess 
as a function of E~, for the same 6-cj) regions. Solid (dashed, dotted, dash-dotted) lines correspond to panels at different <j>, 
starting counterclockwise from <j> = 0° . 
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Binning plot of standard events >100 GeV 
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FIG. 7: Upper panel: Distribution of CTBCORE (the prob- 
ability that the direction estimate is good) and CTBBestEn- 
ergyProb (the probability that the best energy chosen from 
the two energy estimators is correct) for Galactic center line 
events (in red), compared to the distribution of > 100 GeV 
standard events. Lower panel: First tracker layer to show 
evidence of a particle hit for the best track reconstruction. 
Tracker layers are 0-17 where is closest to the calorimeter, 
and 6-17 (2-5) corresponds to FRONT- (BACK-)converting events 
(tracker layers and 1 have no conversion foils). The dark 
(light) gray bars show the average over all > 100 GeV stan- 
dard events ('dirty' sample defined as SOURCE-CLEAN events 
with \b\ > 5°); the red bars show the distribution for the GC 
line. 



distribution is in excellent agreement with the statistical 
expectations (shown by the black dotted line) within the 
±lcr error bars, showing no indication for the presence 
of 'hotspots'. We obtain the same result when reducing 
the size of the regions, shifting them by a common offset, 
or selecting random circular regions all over the sky with 
the requirement that ~ 100 events are included. 



4- Hypothesis: The observed signal is variable. 

An interpretation of the 130 GeV excess in terms 
of dark matter annihilation requires steadiness of 
the source; a strong variability could indicate (time- 
dependent) instrumental effects. In Fig. |H] we plot how 
the TS value of the GC signal evolved over time (assum- 
ing E 7 = 129.8 GeV HH). The dark red line corresponds 



FIG. 8: Histogram of TS values observed in partially over- 
lapping 6° x 6° regions along the Galactic disk, centered on 
6 = 0° and \£\ = 13.5°, 16.5° . . . 178.5°. We show results for 
different line energies from 90 to 170 GeV (dotted lines), as 
well as all TS values combined (black solid line); in the latter 
case, the black dotted line is the theoretical expectation (a 
0.5xi = o + 0.5Xfc=i distribution), the yellow band shows the 
±1ct errors. The tail of the observed distribution looks as 
expected. 



Time evolution of signal 




2010 2011 
time [y] 



FIG. 9: Time evolution of TS values. In dark (light) red 
we show results for the GC region from Tab. |TT] (respectively 
Reg4 from Ref. 21]), in blue we show the evolution of the 
Earth limb line. The gray line indicates for comparison the TS 
value obtained for the Inner Galactic plane, where no signal 
is observed. Jumps in the Earth limb line significance are 
related to times when Z roc k> 52°, cp. Fig. 1161 



to the GC region from Tab. [TT1 the li ght red line shows 
the results for region Reg4 from Ref. 21]. We compare 
them to the time-evolution of the suspicious Earth limb 
line (blue) and the time evolution of the TS value ob- 
tained from the Inner Galactic plane (gray). The 'GC 
region' curve appears to have grown most strongly be- 
tween March 2011 and February 2012, and falling during 
last few months; on the other hand, the signal observed 
in the larger Reg4 does not show the same behaviour and 
appears more steady. In all cases, the open circles indi- 
cate an event between 120 and 138 GeV. The curves show 
no strong sign for a variability of the line features. Note 
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FIG. 10: Incidence angle 9 vs. log E for Z > 110°, with FRONT 
(dark green) and BACK (purple) events. A blue box indicates 
the region 30° < 9 < 45° and 125 GeV < E < 135 GeV, 
where an excess of events appear. The vast majority of limb 
events are at 9 > 60° because the telescope seldom points 
more than 50° from zenith, and the limb events are mostly at 
Z > 110°. 



that the number of GC region line events in the first and 
second half of the observational period are respectively 
12 and 14. 
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III. EARTH LIMB PHOTONS 

In this section, we follow up the weak excess at ~130 
GeV in the limb photons (Fig. [B]) and characterize the 
nature of the apparent excess more precisely. We show 
that the critical incidence angles are not (solely) respon- 
sible for the GC excess. We then consider the possibility 
that an energy mapping error could redistribute events in 
energy, thus making a spectral feature, and explore the 
various parameters of the GC line events and the limb 
bump events. 

In Fig. 1101 we plot incidence angle 9 versus energy 
for limb events (Z > 110°). For 125 GeV < E < 135 
GeV, we find an excess of events with 30° < 9 < 45°, 
as already indicated in Fig. [6] This 9 range contains 
6% of the limb events for E > 100 GeV. The 129 GeV 
bump is clearly visible in the spectrum plotted in the 
bottom panel of Fig. II H and decreases for larger ranges 
of 9 (other panels). The bump appears equally in FRONT 
and BACK converting events. 



A. Energy mapping error: a model for the limb 
bump 

As we discussed above in Sec. II, it is very unlikely for 
the GC line to be caused by extra events (from either 
photons with much higher energies or CR background); 
the same is true for the line feature in the low incidence 
Earth limb events. Given these difficulties, we consider 



FIG. 11: Histogram of Earth limb (Z > 110°) events vs. 
log(iJ) for various ranges of 9. E — 129 GeV is indicated by 
red dotted line. Note the excess in the 30° < 9 < 45° inci- 
dence angle bin, which contains about 6% of the limb events. 
FRONT and BACK converting events are overplotted with dark 
green and purple lines respectively. The 130 GeV excess ap- 
pears equally in FRONT and BACK converting events. 



the possibility that the Earth limb bump results from 
an energy mapping error, and assume that for some un- 
known reason it only affects low incidence Earth limb 
events (and potentially the GC). 6 

We propose a simple model, in which the mapping from 
true energy to reported energy, E(E t ), is linear except for 
a bump near some reference energy. Smooth low-level 
perturbations over large energy scales are not relevant 
here, and could be absorbed in the effective area cali- 
bration. In order to include a small-scale bump in the 
response, we introduce a local compact perturbation in 
the form of a Gaussian. It is convenient to work in log- 
arithmic quantities, so we take x — log Et,y — \ogE, 
and 



Aa exp 



1 (x-X )' 



2a 2 



(2) 



6 Similar to anomalies in the effective area, one would expect that 
such an error should affect all regions of the sky, in contradiction 
to the observations. 
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FIG. 12: Upper panel: function mapping true energy E t to 
reported energy E (see Eqs. [2] and [4]). Lower panel: The 
effect of this mapping on a spectrum dN/dE ~ E~ 2S ', as in 
Eq. [3] (red line) and also for mock data (black histogram). 



where A is a dimensionless amplitude of the bump (— 1 < 
A < 1 is required for monotonicity of y{x)), xq is a ref- 
erence log energy, and a is the width of the bump (see 
Fig. [l~2|) . The effect of the distortion is to change the 
true spectrum dN/dx = dN / dlog(E t ) into an observed 
spectrum 



dN _ dN (dy 
dy dx \ dx 



with 



dx 



Aaexp [ - 



(x — Xq) 2 \ X — Xq 



(3) 



(4) 



Note that that the extreme values oidy/dx — 1±A occur 
at x — xq = ±er and at y — xo(±l — A)a. Assuming the 
true limb spectrum is a power law, we may apply this 
factor to obtain a model spectrum, and maximize the 
Poisson likelihood of observing the data given the model. 

In Fig. [T3l we fit the energy mapping model to the 
Earth limb data for various ranges of inclination angle. 
We find a 4.7a excess of 30° - 45° limb photons at 129 
GeV, with no significant excess at 0° — 30° or 45° — 60°. 



FIG. 13: Fits of the energy mapping model to limb data for 
various ranges of inclination angle 8. The vertical red dotted 
line corresponds to 130 GeV. The test statistic (2AlnL) for 
the best fit model (green line) relative to the null hypothesis 
(red line) is given, along with the significance, expressed in 
"sigma" including a penalty for the 3 additional degrees of 
freedom. The deviation from linearity is only significant in the 
30° < 9 < 45° panel, but not in events with other incidence 
angles. 



B. The Earth limb line and correlations with the 
GC signal 

As shown in the top panel of Fig. [T31 fitting the Earth 
limb events at incidence angles 30° < 9 < 45° with a 
monochromatic line at 129 GeV instead of an energy 
remapping model yields a local significance of only 2.9cr 
(adopting an energy range from 80 to 210 GeV like above 
in Fig. [3|. However, a further tuning of the 6*-range yields 
significances up to 4.1ct (for 25° < 9 < 53°; central panel 
of Fig. IT4j) . but this comes with an additional number of 
trials. In any case, the overall statistical significance for 
a line in the 9 < 60° Earth limb data is above 3er. For 
comparison, the bottom panel of Fig.fMlshows a fit to the 
Galactic center energy spectrum without incidence angles 
30° < 9 < 45°. The GC excess is not removed by this 
cut, which would have indicated a spurious signal. Even 
when removing all events with 25° < 9 < 53° from the 
GC region (from region Reg4 [ll|), we obtain TS = 5.1 
(TS = 10.1) for the Galactic center signal, whereas the 
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FIG. 14: Same as Fig. [3] but for suspicious limb data subset 
with 30° < 9 < 45° (top panel), the 25° < 9 < 53° subset 
tuned to give the largest TS value (central panel), as well as 
for the GC events without the problematic 30° < 9 < 45° 
range for comparison (bottom panel). 



Earth limb line completely disappears. 

The Earth limb line events are distributed all over the 
sky, as expected (Fig. IT5|) . The arrival time of these 
events is concentrated during periods of high rocking an- 
gle, because it is geometrically impossible to see limb 
events at 9 < 45° in normal survey mode (Fig. [16|) . As 
already mentioned in Section UH the distribution of (8, (f>) 
vs. each other and vs. time and longitude are as expected 
(Figs. [T]and |31). In short, none of these tests reveals sus- 
picious trends or correlations, or indicates in any way 
that a systematic error in detector coordinates could map 
events specifically onto the Galactic center. 




FIG. 15: Earth limb events as function of Galactic coordinates 
t and b. The majority of high-incidence limb events appear 
near the orbital pole, which precesses around the celestial 
pole. This pattern is expected from the observing strategy. 
The GC line (red) and Earth limb line (blue) events are shown 
for comparison. The Earth limb line events do not originate 
in the Galactic center. 
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FIG. 16: Time histogram (1-day bins) of GC line (red) and 
Earth limb line (blue) events (upper panel) and all limb events 
(lower panel). The Earth limb line events are only observed 
at high rocking angles that occur during occasional pointed 
observations. The survey mode rocking angle was changed 
from 35° to 50° at 400 days (lower panel). 
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IV. DISCUSSION AND CONCLUSION 

In this paper, we search the publicly available LAT 
data for any trends or correlations that might indicate 
an instrumental origin for the spectral feature at ^130 
GeV towards the GC [H [H Hf. After addressing the 
following concerns, we find no evidence that the 130 GeV 
feature is spurious. 

The Galactic center is bright and has a hard spectrum. 
On degree scales, the GC surface brightness is less than 
a factor of 2 brighter than the inner Galactic plane. The 
inner plane provides an order of magnitude more pho- 
tons, and shows no sign of a 130 GeV bump. Even larger 
samples (all limb photons, all non-limb non-GC photons) 
also show no significant signal. The GC has a hard spec- 
trum, and at TeV energies the GC is brighter than the 
surrounding plane, but even if all events above 300 GeV 
(assuming a hard spectrum dN/dE ~ E~ 2 ) were incor- 
rectly remapped to 130 GeV, there would not be enough 
photons to explain the GC signal. 

Observations of the GC have a restricted range of in- 
strumental incidence angles. It is true that the survey 
strategy, orbital precession, and solar panel alignment 
cause a non-trivial mapping of GC events onto (9, <fi) as a 
function of time of year. Specifically, they occur near the 
x — z plane (4> — or 180) when the Sun is near the GC 
or anti-center. However, the GC line events are drawn 
from the full range of 9, 4>, and t. The limb line events 
are broadly distributed in 9, 0, and £, with times corre- 
sponding to pointed observations at large zenith angle. 
There is no evidence that the distribution of line events 
deviates from expectations. 

There are excess line events in the limb data for some 
incidence angles. For a small subset of the limb data 
with large rocking angle (when the limb may be seen at 
small incidence angles) and a particular incidence angle 
range around 30° to 45°, we find a marginally significant 
130 GeV feature (above 3a). The majority of events 
with incidence angle 30° to 45° are not from the Earth 
limb, and we find no 130 GeV feature in this much larger 
sample of events at these incidence angles. If the limb line 
events are an artifact, they must conspire to only appear 
when the LAT is positioned at high rocking angle. 

The bump in the limb data might result from an energy 
mapping error. We propose a simple model for an error in 
the mapping from true photon energy to reported energy. 
This model reproduces the shape of the limb line feature 
at 130 GeV and the dip at slightly higher energy, and 
has a local significance of 4.7a. A modest amount of 
additional limb data would tell us if the limb feature is 
a statistical fluke. If the limb feature persists, it raises 
serious concerns about the Pass 7 processing of E > 100 
GeV events. 

Additional limb data are available from the commis- 
sioning period. The Launch & Early Operations (LEO) 
data were taken during the first 60 days of the mission. 
Combined with a dedicated Earth-limb observation in 



September 2008, this provides ^250 hours total livetime 
on the Earth limb [491 ]. 

With Pass 6 diffuse class events, 49] has analyzed 
the spatial morphology and the energy spectrum of the 
Earth limb sample, which contains 218 photons above 
100 GeV and 16 photons above 500 GeV. The energy 
spectrum is a power-law with spectral index 2.79 ± 0.06 
for 3-500 GeV photons (Fig. 2 of |49|), which is consistent 
with the primary cosmic-ray spectral index 2.75 ± 0.03. 7 

The spectrum of the Earth limb photons provided by 
[49l | does not show any significant feature at 130 GeV. If 
improved processing of the limb photons does not repli- 
cate the line or the " energy mapping error" we found for 
a subsample of the limb photons during the normal sur- 
vey mode, it can be dismissed as a statistical fluke. If 
it reappears, a deeper investigation into its cause will be 
necessary. 

Even then, it is a challenge to understand how such an 
instrumental feature could be mapped so precisely onto 
a localized region within 5 — 10° of the GC. The GC 
is not near the path of the orbital pole, nor its axis of 
precession. The orbital phase, precession, Earth's orbit, 
and time of year are all well mixed by the few x 10 4 orbits 
and 25 precession cycles over 1500 days. We have shown 
that the events in question are drawn from every part 
of event and spacecraft parameter space available in the 
public files. 

In summary, we find no significant instrumental sys- 
tematics that could plausibly explain the excess Galactic 
center emission observed at 130 GeV. 

Note added: During the final stages of this work we 
became aware of another group discussing instrumental 
indications in the Earth limb data [55[. 
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7 The Earth limb photon above 10 GeV from CR interactions in 
the upper atmospheric layers do not suffer large energy loses and 
the cosmic-ray primaries at this energy are unaffected by the 
Earth's magnetic fields. The Earth limb photons should have a 
spectral index close to that of the primary cosmic rays [49| . 
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FIG. 17: Distribution of events in Earth longitude and lati- 
tude (upper panel) and satellite zenith RA,dec (lower panel). 
In this and following figures, Galactic center line photons 
(red), Earth limb line photons (blue) and Earth limb photons 
for all incidence angles (black) are shown. The avoidance of 
the SAA leaves a hole in the upper panel. 



FIG. 18: RA and Dec of spacecraft Z axis (boresight direc- 
tion) and X axis (Solar panel direction). 



http : //idlastro . gsf c . nasa. gov). 



V. APPENDIX 

In this Appendix, we compare the distribution of the 
Galactic center line events with the Earth limb photons 
in various projections of the event parameter space, and 
search for any unexpected behavior. As above, red points 
represent the Galactic center line photons and the blue 
points represent Earth limb line events. Earth limb pho- 
tons with 135 GeV > E > 125 GeV at all incidence angles 
are shown in black. None of the Figures (Figs. rT?ll2"2")) 
show any sign of unexpected behavior. 
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FIG. 19: Zenith angle vs. Earth azimuth angle for limb pho- 
tons. As expected, the 8 > 60° limb photons observed in 
survey mode are seen predominantly to the north and south 
azimuth directions, i.e approximately perpendicular to the or- 
bit direction. The blue points are the Earth limb line events. 
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FIG. 22: All events as function of zenith angle Z and incidence 
angle 6. The change of the rocking angle from 35° to 50° in 
2009 is clearly visible. None of the GC line events is observed 
near the horizon. 



FIG. 20: Events as a function of orbital precession phase 
(upper panel) and time of year (lower panel). Note that there 
are only about 6 months of the year during which the 38 Earth 
limb line events (blue) are seen. 
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FIG. 21: The rocking angle of the Galactic center line events 
(positive values indicate a rock angle toward the north from 
zenith. 
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